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O Abstract 

<N We measure the quantum fluctuations of a pumped nonlinear resonator, using a superconducting 

^ artificial atom as an in-situ probe. The qubit excitation spectrum gives access to the frequency and 

^ temperature of the intracavity field fluctuations. These are found to be in agreement with theoretical 

predictions; in particular we experimentally observe the phenomenon of quantum heating. 

PACS numbers: Valid PACS appear here 
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A resonator in which a nonhnear medium is inserted has rich dynamics when it is driven 
by an external pump field [1]. In the case of a Kerr medium, the field amplitude inside 
such a Kerr nonlinear resonator (KNR) switches from a low to a high value when the 
pump reaches a certain threshold, a phenomenon known as bistability in optics [2j and 
bifurcation in the microwave domain [3j. This hysteretic transition between two dynamical 
states is a stochastic process triggered by fiuctuations around the steady-state pumped 
oscillations. Given the potential applications for low-power all-optical logical elements such 
as switches and transistors [4j, a quantitative understanding of these fiuctuations governing 
the sharpness of the transition, and therefore the device performance, is highly desirable 
[HI E]. In a KNR operated in the quantum regime [7j, field fiuctuations are mainly due to 
spontaneous parametric down-conversion (SPDC) [814T0] of pairs of photons at the pump 
frequency Up into pairs of photons, one at the characteristic pumped-KNR frequency Uc and 
the other at the complementary idler frequency = 2ujp — Uc- The mode at therefore 
acquires a thermal population with an effective temperature that depends on the pump 
frequency and power even when the electromagnetic bath to which it is coupled is at zero 
temperature, a central theoretical prediction known as quantum heating that remains to be 
tested [ZmiHIS]. 

Up to now, experiments in the optical or microwave domain have only measured the 
spectrum of the field radiated by a pumped KNR jT3]. In this work, we access the intra- 
resonator field fiuctuations by inserting a two-level system (TLS) inside the KNR and using 
it as an absolute spectrometer and thermometer, as illustrated in Fig. [T^. The experiment 
is performed in the microwave domain at millikelvin temperatures using superconducting 
Josephson circuits. The KNR is a coplanar waveguide resonator with an embedded Joseph- 
son junction [T^ FIB], and the TLS is a transmon qubit [T7] (see Fig. [Tjo). By measuring the 
qubit absorption spectrum while pumping the KNR, we obtain a quantitative agreement 
with theoretical predictions on quantum heating. 

We start with a rapid summary of the relevant theoretical results regarding the quantum 
fiuctuations of a KNR - more details can be found in |5l [71 [TT | [T3 | [rSVEO] . The KNR, driven by 
a field of power Pp and frequency Up, is modelled in the frame rotating at Up by the Hamilto- 
nian 

Hc/h = Apa^a + ^a)'^a? + ^a)^a? + [itpo) + 'i{.c). with Ap = uJc — (jJp, uic the KNR resonance 
frequency in its linear regime, d and the KNR field annihilation and creation operators, 
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ep = a/ /t-Pp/ hUp the driving amplitude, k, the resonator energy damping rate, and K and K' 
the Kerr nonhnear constants derived from circuit parameters [T61I21) . The steady-state solu- 
tion for the dimensionless cavity field amplitude a is obtained from the corresponding master 
equation yielding (Q^ + Kla]"^ + K'\a\* — i^) a = —ie^ with VL = 2QAp/ijJc the reduced 
pump frequency [T8| [20] . and Q = Uc/ n the resonator quality factor. For Q > y/S and Pp 
larger than the critical power P^. = hwp this equation admits two stable solutions aL,H 
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FIG. 1. (a) Optical analogue of the experiment: a Fabry-Perot cavity with resonance frequency Uc, 
quahty factor Q, and Kerr non-linearity is pumped at frequency ujp. Photons are produced at 
the cavity dressed mode frequency Cjc and at the idler frequency by SPDC of pump photon pairs. 
An in-situ atomic probe is used as an absolute spectrometer and thermometer to characterize 
this quantum heating of the dressed mode, (b) Actual experiment: the cavity is a microwave 
superconducting A/2 coplanar waveguide resonator, the Kerr medium a Josephson junction, and 
the probe a transmon artificial atom, (c) Stability diagram of the resonator in the pump power 
Pp and frequency Up plane (in reduced units, see text) : L, i7, and B denotes regions where the 
amplitude of the driven oscillations is low, high, and bistable, respectively. P4. is the pump power 
at which the field switches from L to H when Pp is increased, (d) Dressed mode frequency uJq when 
increasing Pp at = — 3 (along the arrows of panel c). 
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corresponding to metastable dynamical states L and H of respectively low- and high- ampli- 
tude [18]. In this bistable and hysteretic regime, the transition from L to H occurs abruptly 
when ramping up the pump power at the bifurcation threshold P^{Q). The corresponding 
stability diagram is shown in Fig. [T]:. 

We consider here the quantum fluctuations of the intra-resonator field around its steady- 
state value a, in the regime where they are too weak to induce switching to the other 
dynamical state within the experiment duration. This justifies linearizing He around a, by 
writing a = a + a and keeping only terms quadratic in a. Following j2Q], this linearized 
Hamiltonian can be diagonalized by introducing a new operator a = ^a + va\ and rewritten 
as Hi = Apd^a with Ap = sign{A)^/B, A = + 2K\a\'^ + 3K'\a\^, and B = A"^ - {K + 
2K'\a\'^)^\a\'^. In the laboratory frame, intra-cavity field fiuctuations are thus described as 
excitations of a harmonic dressed mode of resonance frequency a)c = Up + Ap that depends on 
the pump amplitude and frequency |22| . We note the eigenstates of this effective oscillator 
\n) = (a"!")" |o) / Vn). . As shown in Fig. [iji for the case where K, K' < and Ap > as in 
our experiment, the dressed frequency cDc is equal to Uc at low pump power, then decreases 
when Pp is increased, reaching Up when Pp = P+, which causes the field to jump to its 
high amplitude value and uj^ correspondingly to jump well below Up. The dressed mode 
is damped at the same rate k, as the KNR, but towards an equilibrium steady-state at a 
finite effective temperature T^s corresponding to a mean number of excitations (n) equal to 
|z/p, even if the bath physical temperature is zero |20]. Physically this thermal population 
is caused by SPDC of pairs of pump photons at Up into correlated photons at frequencies 
cDc and Ui = 2ujp — Uc, emitted in the dressed mode and in the measuring line respectively; 
the apparent thermal character of the intraresonator field is obtained when neglecting the 
correlations between the 0^ and the Ui photons. Note that this analysis is only valid if 
-B > 0, a condition verified sufficiently far from the bifurcation threshold, as is the case here. 

In our experimental test of these predictions, the KNR is a superconducting coplanar 
waveguide resonator including a Josephson junction with frequency Uc/27i = 6.4535 GHz, 
Kerr constants K/2tt = —625 kHz and K' /2tt = —1.25 kHz, and damping rate k,/2tt = 
10 MHz (see Fig. [l|l) [16j. It is capacitively coupled with strength g/2-7T = 44 MHz to a 
superconducting qubit of the transmon type with frequency Uge/'2iT = 5.718 GHz. Because 
of the large qubit-resonator detuning, their interaction can be described in the so-called 
dispersive limit p3]. Resonator and qubit can be driven by microwave pulses applied to the 
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FIG. 2. Measurement of the field radiated by the KNR pumped at = 2.8. (a) Experimental setup: 
the pump tone, sent through a line with attenuators (A), is switched on and off every 2.5ms; the 
reflected and radiated signals are amplified and mixed down at Wp, yielding two quadratures I{t) and 
Q{t) that are digitized and Fourier transformed; the total voltage spectral densities 5on,off(^'^ = 
|c<j — (^pl) for the two pump states are then computed, (b) Radiated noise spectrum Son — •S'oFF 
for Pp = — llOdBm. (c) Radiated noise spectrum as a function of Pp. The dashed curve is the 



calculated Ar 



a;p(Pp). The vertical dotted line indicates the cut shown on panel b. 



resonator input. The qubit can be readout in a single shot by driving the resonator close to 
its bifurcation threshold P+{^) in order to map the qubit states \g) and |e) to the L and H 
states, respectively |24| . 

A first method to investigate the field fluctuations of the KNR pumped at uj^ is to measure 
the field radiated into the measurement line at frequencies d)c(Pp) and 2ujp — udPp), as 
reported recently [HI Its noise spectrum, measured using the setup of Fig. ^ \2Q\ . 

shows a peak at a Pp-dependent frequency (see Figs. |2|3-c). All the sample parameters 
being known from earlier measurements [IE], the theoretical a)c(Pp) curve can be computed 
without adjustable parameters and is shown in Fig. |2]3. Except for avoided crossings at 
Ap/27r ^ 40 and 70 MHz of unknown origin, the agreement is quantitative. This shows that 
the additional noise is indeed generated by the pumped KNR at ujc{Pp). However these 
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FIG. 3. Sideband spectroscopy of the qubit while pumping the resonator with varying power Pp at 
0, = 3.9. (a) Experimental setup: a spectrocopic pulse is applied at coq with power Pq during the 
pumping; the pump and spectrocopy tones are then switched off and the qubit is measured using 
a third pulse at cji- with a power close to the bifurcation threshold P+. The measured phase of 
the reflected pulse yields the qubit state and repeating the sequence yields the qubit excited state 
probability p|g^.(b) Measured (top) and calculated (bottom) 2D plots showing the p|g^((jjq) spectra 
as a function of Pp (bottom axis), also expressed in units of the bifurcation threshold P-|_ (top 
axis), (c) Experimental (dots) and analytical (lines) spectra at a power Pp = — 96dBm (dashed 
lines in b), for powers Pq = -131, -120 and -115 dBm (from left to right), (d) Energy diagrams 
showing the transitions involved in the red sideband peak (Stokes), central peak, and blue sideband 
(anti-Stokes). 



measurements are unable to determine the effective temperature of the mode uJc from which 
the measured photons are leaking, which is the key quantity of the theory discussed above. 

A two-level system such as the transmon qubit is an ideal in-situ probe of quantum 
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heating, because it acts as an absolute spectrometer and thermometer for a given quantum 
noise source |27J. We measure the temperature of the a)c mode while keeping the qubit 
and KNR far detuned by using a method called sideband spectroscopy, which is routinely 
used in ion-trapping experiments [29| [30] and has also been applied recently to mechanical 
oscillators [31\. Indeed, starting from the system in state \g,h) it is possible to drive a 
e,n + 1) by irradiating the qubit at frequency ujge + Coc (so-called Stokes 



transition to 

sideband), or to e,n — 1^ provided n > by irradiating the qubit at \ujgc — (so-called 
anti-Stokes sideband) [32VI34) . In the case of a transmon qubit, these transitions need to 
be driven with two photons of arbitrary frequency provided their sum (resp. difference) 
satisfies the Stokes (resp. anti-Stokes) sideband resonance condition ^35j. Given the sideband 
transition matrix element dependence on n, one can show that the average photon number 
(n) = l/[exp{hu}p/kTes) — 1] is equal to r/(l — r) with r G [0,1] the anti-Stokes/ Stokes 
sideband peak amplitude ratio |30]. Measurement of the qubit absorption spectrum, yielding 
r, therefore corresponds to a direct and absolute measurement of Tgff. 

Sideband spectroscopy is performed using the setup shown in Fig. |3^. Once a steady-state 
pump field is established in the resonator, a spectroscopy pulse is applied to the qubit at 
fixed power Pq and varying frequency Uq. The pump tone provides one of the two photons 
needed to drive the sideband transitions, and the spectroscopy tone provides the second 
photon whenever Uq matches the Stokes (resp. anti-Stokes) sideband resonance condition 
Wq + Wp = Uge + 0Jc{Pp) (resp. Up — Uq = ojc{Pp) — Wgo)- The experimental sequence ends 
by reading out the qubit state 200 ns after both pulses are switched off, long enough for the 
KNR field to decay but shorter than the qubit relaxation time Ti ^ 700 ns [16j; repeating 
this sequence ^ 10^ times yields the qubit excited state probability p\f,y 

Typical data are shown in Fig. [3]:3 at = 3.9 and Pp = —96 dBm such that the KNR is in 
the high oscillation amplitude state H. At low spectroscopy power Pq, only one Lorentzian 
peak is visible, corresponding to the qubit frequency ac-Stark shifted by the steady-state 
intraresonator field with mean photon number (nn) = Itt/zP at Up [TBI 136] . Increasing the 
spectroscopy power, we observe the appearance of two satellite peaks around Uge, with a 
separation of 31 MHz that closely matches the value of Ap = Up — ujdPp), already known 
without any adjustable parameters as explained above. When Pp is varied, this separation 
also quantitatively varies as expected from the Uc dependence on Pp, as shown in Fig. |3|3. 
This establishes that the satellite peaks are indeed the sideband transitions. The anti-Stokes 
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FIG. 4. Quantum heating of the pumped resonator mode Cjc above bifurcation: Quantity — r) 
with r being the ratio of blue to red sideband peak heights deduced from the experimental (squares) 
and theoretical (dots) qubit spectra of Fig. [sj), as a function of the pump power Pp (bottom axis). 
The calculated average number (nn) of photons in the pump mode is indicated on the top axis. 
In a simple theoretical picture r/(l — r) is the average photon number (n) in the Cjc mode, shown 
here as a dashed and a solid line for the L and H resonator states. The corresponding effective 
temperature Tgfj is indicated on the right axis. 




sideband being observable and of smaller amplitude than the Stokes sideband indicates that 
the temperature of the dressed mode is finite, as discussed in more details below. 

To be more quantitative, we have performed a detailed theoretical analysis [20] of the 
coupled qubit-KNR system that will be presented elsewhere [37], and which yields analytical 
approximate expressions for the qubit sideband spectrum. The predictions, calculated with 
a global attenuation factor on the spectroscopy power Pq as the only adjustable parameter, 
are also shown in Fig. |3}3-c for different Pq; they agree quantitatively with the data. Since 
these calculations are done at zero bath temperature, this is a first clear indication that the 
population of the mode at a)c is only due to SPDC and is therefore of quantum origin. 

Assuming that the field in the dressed mode uj^. has the statistics of a thermal field as 
predicted theoretically, we now extract an experimental occupation number of the dressed 
mode from the relative height of the two sideband peaks, and translate it into an effective 
temperature Tofj. For this, each spectrum of Fig. [Sjo is fitted to a sum of three Lorentzians 
of adjustable frequency, width and height, yielding the ratios r of the anti-Stokes to Stokes 
sideband height. Figure |4] shows the comparison between the experimental occupation num- 
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ber r/(l— r) and (n) calculated without any adjustable parameter. The data agree fairly well 
with the prediction, demonstrating that the average photon number in the dressed mode is 
indeed |z/p as predicted by theory (see chapter 7 in [Ij). More precisely, we estimate that the 
residual thermal field at Cjc cannot be responsible for more than 20% of the observed signal. 
Our data show that the dressed mode occupation number and temperature are maximum 
near the bifurcation threshold Pj^ and decreases with increasing pumping power and occu- 
pation number (wh) at uj^. This constitutes additional evidence that the measured thermal 
field is not due to a trivial heating caused by the microwave pulses. We also stress that these 
results do not rely on any calibration of the measurement lines or temperature. Finally, we 
also show in Fig. |4] the quantity r/(r — 1) derived from the analytical expression that yielded 
the theoretical spectra shown in Fig. [3] it is in agreement with both the experimental data 
and the simple formula for the KNR effective temperature. Additional data can be found 
in the Supplementary Information [20j taken on another sample with different parameters, 
with which similar results have been obtained. 

In conclusion, we have probed the quantum fluctuations of a pumped nonlinear resonator 
with an embedded superconducting qubit, bringing experimental evidence of quantum heat- 
ing. Future directions include establishing the link between quantum heating and the switch- 
ing rates at the bistability threshold [llj, and testing the thermal character of the dressed 
resonator mode by performing its quantum state tomography with the qubit [38j. In general 
our experiments demonstrate that detailed and quantitative tests of all theoretical predic- 
tions regarding nonlinear resonators in the quantum regime are enabled by the progress of 
superconducting circuits. 

We acknowledge useful discussions with M. Dykman and within the Quantronics group, 
as well as technical support from P. Senat, P.-F. Orflla, T. David, and J.-C. Tack. We 
acknowledge support from NSERC, the Alfred P. Sloan Foundation, CIFAR, the European 
project SOLID, the C'Nano project QUANTROCRYO, and the Australian Research Coun- 
cil via the Centre of Excellence in Engineered Quantum Systems (EQuS), project number 
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SUPPLEMENTARY INFORMATION FOR "QUANTUM HEATING OF A NON- 
LINEAR RESONATOR PROBED BY A SUPERCONDUCTING QUBIT". 

THEORY OF QUANTUM FLUCTUATIONS OF A NONLINEAR RESONATOR 

Our goal here is to give a self-contained account of the theory of quantum fluctuations of 
a nonlinear resonator. 



Driven Kerr shifted cavity: linearized description 

Consider a Kerr nonlinear cavity with frequency Uc, Kerr coefficients K and K', and 
driven by a field at frequency oUp and amplitude e^. The Hamiltonian describing the driven 
cavity is 

K K' 

He = huJcO^a + h—a^'^a^ + h—a'^a^ + 2hepSm{ujpt){a + a^) (1) 

with a and the KNR field annihilation and creation operators. 

In a frame rotating at uip and under the rotating wave approximation, we find 




where we have defined Ap — Uc — ojp. Allowing the drive amplitude to be complex in this 
equation allows for different choices of drive phase. Including cavity damping k — uic/Q {Q 
being the cavity quality factor), the master equation in this frame is then 

-i[Hc,p\/h^ KD[a\p, (3) 

where 

V[c]p = cpc^ - c^cp/2 - pc^c/2. (4) 

In the limit in which the linearization is appropriate, the cavity is not too far from being 
in a coherent state. It is then useful to define new operators a, for the cavity such that 

a — a-\-a, (5) 

where o; is a complex number which will be shown later to correspond to the coherent state 
amphtude in the resonator. The operator a describes fluctuations about this coherent state, 
and the contribution to the photon number from this state should be small compared to lap 
for the linearization to be a valid approximation. 
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Making this replacement in He, we find the hnearized Hamiltonian 

Hi = h{A + 2K\a\^ + 3A^'|«|^) a^a + h (^^a*^ + K'\a\^a*^^ a^ + h(^ja^ + K'\a\^a^^ 

(6) 

where the value of a is chosen such that all terms linear in a, including those arising from the 
dissipation, cancel (see below). We have also disregarded terms that involve higher powers 
of a and on the grounds that the effects of these terms will be small since they involve 
lower powers of a which is assumed to be large compared to one. 

Bistable steady state amplitudes 

The value of a is found by requiring that, after substituting a = a + a (Eq. |5]) in the 
master equation Eq|3j the terms proportional to a have a zero coefficient. We first give the 
term linear in a originating from the Hamiltonian part of the master equation: 

h{-itl + Aa* + K\a\^a + K'\a\'^a) . (7) 

In the dissipative term of the master equation, the substitution Eq. |5] leads to 

P[a]p = V[a]p — i[ih—{Q*a — aa)), p\/h, (8) 

the last term of which contributes to the pref actor of a. Gathering all the terms linear in a 
leads to a verifying 

iel = Aa* + K\a\^a* + K'\a\^a* + %^a* . (9) 
Introducing the reduced variable f2 = IQA^juc this can be rewritten 

- iep = (n'^ + + K'\a\^ - a, (10) 

as given in the main text. 



Diagonalizing the linearized Hamiltonian 

In order to understand the linearized spectrum of the cavity we diagonalize the linearized 
Hamiltonian 
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Hi^h{A + 2K\a\^ + 3K'\a\^) a^a+h + i^Vl V^^ a^+h (^^a^ + i^Vl 

(11) 

We wish to diagonalize this by means of a Bogohubov style squeezing transformation by 
defining a new "lowering" operator 

a = /xa + ua^ (12) 



and its corresponding conjugate 
Simple calculation shows that 
if and only if the condition 
holds. 

Consider now the Hamiltonian 



at = i_i*a^ + i/*a. (13) 



[a,at] = l (14) 



;^|2_|ix|2^1 (15) 



Hi = HApo'd. (16) 



We will show that for a suitable choice of /j, and u this Hamiltonian corresponds to Hi plus 
a constant. 

Without loss of generahty we will be able to choose // to be real. This constraint will hold 
identically if we define /i and u in terms of hyperbolic functions /i — coshr, u — e^*^sinhr. 
Again without loss of generality, we may also require that r > and hence sinhr > 0. 

The most straightforward way to proceed is to write the new Hamiltonian in terms of the 
original operator and match coefficients : 

Hi = hAp [ii*a^ + v*a) [jia + va)) (17) 
= Mp + a^a + Mp (/xVa^^ + /xi/*a^) + hA^\v\^. (18) 

We can see that this is equal to Hi up to a constant as long as we can choose 

Ap cosh(2r) = A + 2K\a\^ + 3fs:V|^ (19) 
2_,iinh(2r)e-2'^=(| 



^-ApSinh(2r)e-2^^ = + ^Vl' ) (20) 
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Notice that the first equation determines the sign of Ap since cosh is positive: Ap must 
have the same sign as A + 2ii'|ci;|^ + SiC'lal^. 

Looking at the second equation, we will define ol = |Q;|e*^= and then we can conclude 
that exp(— 2i^) = ±exp(— 2i^c). Recalling that we have chosen r > implying that sinhr 
is positive, the plus sign is obtained when Ap has the same sign as K and the minus sign 
otherwise. In the former case we have 9 — 9c and in the latter 9 — 9c + n/2. In either case 
the squeezing axis rotates with the phase of the cavity amplitude. 

By taking the ratio of the absolute values of these equations we find an expression for r 



tanh 2r = 





\K + 2K'\ 


\a\ 


21 


\a\ 


2 


|A + 2ii'| 


\a\ 


2 


+ 3K' 


\a\ 


4 



(21) 



while the magnitude of Ap can be inferred from squaring and substr acting the two equations: 

A^ = (A + 2K\a\^ + 3K'\a\^y - {K + 2K'\a\^f\a\\ (22) 

The regime of interest for us is where (A + 2ii'|Q;p + 3K'\a\^f - {K + 2K'\a\'^)^\a\'^ > 
in which case Ap is real and the dynamics of the KNR fluctuations can indeed be mapped 
onto those of a damped harmonic oscillator of frequency Ap in the rotating frame. This is 
the case in most of the KNR phase space, except very close to the bistability threshold, in 
which case the fluctuations are instead amplifled and the system behaves as a parametric 
amplifler of large gain. In this work we only focus on the regime where Ap is real. 

In this regime we can thus interpret the system as a cavity of frequency Uc with = 

Ap + LUp. 

Approximate master equation in the linesirized description : quantum heating 

In the regime where k <^ (Dc, we can move into a new interaction picture at frequency 
ujc- In this interaction picture, the Hamiltonian is zero and only the dissipative terms in the 
master equation contribute. We have so far written these terms only as a function of a but, 
when we express them in terms of d using the inverse transform mentioned above and move 
into the interaction picture, we flnd that some terms have a time dependence at frequency 
2Ap. In the limit Ap ^ k in which we are working (corresponding to the resolved-sideband 
limit), one can make a rotating wave approximation and average these terms to zero. The 
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resulting master equation is 



p = fi;(|zy|2 + l)V[a]p + K\i^\'^V[a'<]p. 



(23) 



This is the master equation for an oscillator coupled to a thermal bath with energy decay 
rate n and thermal expectation number The steady state of the cavity mode is thus a 
thermal state with (d'^a) = |z/p. It may seem strange to have a thermal state as a steady 
state. Physically the quantum state in this interaction picture approaches a squeezed state 
with some periodic time dependence of the phase of the squeezing. This variation is so rapid 
that, as long as we respect the regime of validity of the rotating wave approximation, we can 
replace it with a time averaged state which is thermal. This phenomenon has been coined 
'"quantum heating'" jlHS]. 

ANALYTICAL FORMULA FOR THE SIDEBAND SPECTRUM 

The analytical results presented in Fig. 3 of the paper were obtained using a more com- 
plete theory than that presented above. This theory is based on the one developped in 
Ref. [4J, and will be developped further in Ref. [H]. It includes effects from the resonator 
on the qubit such as Purcell decay, measurement-induced dephasing j6] and dressed dephas- 
ing j7] and decay, as well as Lamb and ac-Stark shift of the qubit transition frequencies. It 
also includes effects from the qubit on the resonator, such as a qubit-state dependent pull 
of the resonator's frequency leading to different pointer states of the resonator for each 
qubit state The theory starts from a master equation description of a qubit coupled 
to a nonlinear resonator driven both with a pump ep and spectroscopy drives. Through 
successive transformations |5], we transform the master equation into a form in which the 
resonator can be adiabatically eliminated, yielding a reduced qubit master equation given 



Here, as is the intra-resonator field created by the spectroscopy drive of amplitude at 
i^s = ^^1,0 ~ close to the qubit transition frequency. The dephasing rate 7^ includes the 



by 



Pq = -AH, Pq] + liV[aJiPq + ^^V[a+]pq + ^P[ajpg, 



(24) 



where 




(25) 
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intrinsic qubit dephasing as well as the measurement-induced dephasing [5|. The up and 
down rates are given by 



(26a) 
(26b) 



where 7|"|. include Purcell relaxation and dressed dephasing and are obtained in Ref. |1] and 
L{5) is a Lorentzian of full-width at half-maximum n given by 

k/2 



m 



/€V4 + (Ap + 5)^ 



(27) 



In this equantion Ap is given by Eq. 22 and /(oj, jj) is a unit-less function of the squeezing 
coefficients and of the distinguishability of the pointer states ai — ao- 



The stationnary solution of the reduced master equation is analytical and yields 



(ni,i),q (71 + n + 272 |^?o«.,o| /(7t + li) 



where 



[(7I + 472 |^oa.,o|^ /(7t + li)) + 



7t 



' cq 



7t + li 



(28) 



(29) 



and 72 = % + (7t + 74.)/2. Equation (28) was used to plot the analytical data presented in 
Fig. 3 of the paper. 

In the low spectroscopy power regime, the effective rates 7t,4.,2 are independant of the 
spectroscopy frequency, and the spectrum yields a single line, centered at 5 = 0, with a 
minimal width of 72. In the high spectroscopy power regime, Us becomes significant and 
the central line is power-broadened. Moreover, the rates acquire a frequency-dependent 
structure. If we focus on the well-resolved sidebands limit, that is 5 > 72, the dominant 
term yields -P(|l)) ~ (ni i)^^. One can then easily compute the blue to red sidebands ratio 
of amplitude simply with 



(ni,i) 



eq 



<5=Ab 



(ni,i) 



eq 



7tlg=A^ 
7tl5=-A 



L(-Ap) + L(A, 



L{-Ap) + L{Ap) \u\^ + L{-Ap) 





V 


2 




2 


+ 1 



(30) 



In the first approximation, we assumed that the up rate 7^" is negligible, while in the second 
approximation, we assumed that L(— Ap) ^ L{Ap), which should be the case in the well 
resolved sidebands limit. 
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DATA FROM ANOTHER SAMPLE 



To complement the data shown in Fig. 4 of the main text, we now present data showing 
sideband spectroscopy of a similar sample, although with different values of the parameters. 
For this sample B, = 6.469 GHz, Q = 1040, and K/2n = 1 MHz. K is determined as 

explained in [6| by fitting the frequency shift (AC Stark shift) of the qubit while the KNR 
is pumped. The pump frequency was chosen to be Up/27T = 6.427 GHz, corresponding to a 
significantly larger Q = 13.5 than in the main text. 

Sideband spectroscopy was performed similarly to what is reported in the main text. 
Data are shown in FigjSj around the bistability threshold. Above the bistability threshold 
the data are similar to the data presented in Fig.3 of the main text. Note the additional 
presence of a faint fourth line at Uge + 2(a)c(Pp) — Wp), that corresponds to a transition from 
\g,h) to e,n + 2^, and of a narrow peak that corresponds to the two-photon excitation 
of the qubit from state l^f) to its second excited state, AG-Stark shifted by the field inside 
the KNR. Below threshold (resonator in state L), the Stokes sideband is also well resolved 
and visible. This was not the case in the sample discussed in the main text because the 
pump frequency was much closer to the KNR frequency in these data, resulting in a smaller 
separation between the Stokes sideband below trheshold and the main qubit peak (since 
well below the threshold Ap ^ Uc — ujp), which made the two lines indistinguishable. Quite 
remarkably, even though the Stokes sideband is clearly visible below the threshold, the anti- 
Stokes sideband is invisible within our detection limit (see Figjs]). Since the ratio between 
Stokes and anti-Stokes peaks directly yields the mean photon number (n), this sets an upper 
bound on the mean photon number (h) present in the mode at Uc which is shown in Figjs] 
for various pump powers (the dependence on pump power is simply due to the fact that the 
height of the Stokes peak also diminishes at small pump power, whereas the noise in the 
data is constant). In particular, this allows us to establish that any residual photon numbers 
that might be present at cDc due to improper thermalization or filtering is below the lower 
point of the dashed blue line in FigjS} namely < 0.04. This is again strong evidence that the 
much larger thermal photon numbers observed above threshold are genuinely of quantum 
origin. 

Another noteworthy feature is that the Stokes and anti-Stokes sideband position with 
respect to the qubit frequency Ug are inverted at the bistabihty threshold, as expected 
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FIG. 5. (a) (color plot) Sideband spectroscopy of sample B on both sides of the bistability threshold, 
i.e. for the resonator in the L state (Pp < — lOSdBm) and in the H state (Pp > — 105dBm). The 
expected frequency of the AC Stark-shifted main qubit peak as well as the Stokes and anti-Stokes 
sidebands are shown as dashed white and black lines, (b) (red dots) Measured r/(l — r), with 
r being the ratio of the anti-Stokes to Stokes sideband peak height. Data points are shown only 
above the threshold, because no measurable anti-Stokes peak is found in our data below threshold. 
We show instead our detection limit, defined as n//(l — ra/), with nf the ratio of the standard 
deviation in p|g^ to the Stokes sideband peak height, above (resp. below) threshold, as a blue (resp. 
red) dashed line. The detection limit depends on the pump power because the Stokes sideband 
height does. The value of (n) calculated without adjustable parameter is shown as full line above 
(red) and below (blue) the threshold. 

from the sign change of Ap (see Fig.l in the main text). Finally, the effective temperature 
extracted from FigjS^ is shown in Fig|5|D and is again in fair agreement with the prediction 
|z/p calculated without adjustable parameter, both below and above the threshold. 
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